Pairwise entanglement and readout of atomic-ensemble and optical wave-packet 
modes in traveling-wave Raman interactions 
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We analyze quantum entanglement of Stokes light and atòmic electrònic polarization excited 
during single-pass, linear-regime, stimulated Raman scattering in terms of optical wave-packet modes 
and atomic-ensemble spatial modes. The output of this process is confirmed to be decomposable 
into múltiple discrete, bosonic mode pairs, each pair undergoing independent evolution into a two- 
mode squeezed state. For this we extend the Bloch-Messiah reduction theorem, previously known 
for discrete linear systems (S. L. Braunstein, Phys. Rev. A, vol. 71, p. 055801, 2005). We present 
typical mode functions in the case of one-dimensional scattering in an atòmic vapor. We find that in 
the absence of dispersion, one mode pair dominates the process, leading to a simple interpretation 
of entanglement in this continuous-variable system. However, many mode pairs are excited in the 
presence of dispersion-induced temporal walkoff of the Stokes, as witnessed by the photon-count 
statistics. We also consider the readout of the stored atòmic polarization using the anti-Stokes 
scattering process. We prové that the readout process can also be decomposed into múltiple mode 
pairs, each pair undergoing independent evolution analogous to a beam-splitter transformation. We 
show that this process can have unit efficiency under realistic experimental conditions. The shape 
of the output light wave packet can be predicted. In case of unit readout efficiency it contains only 
excitations originating from a specified atòmic excitation mode. 

PACS numbers: 42.65.Dr,42.50.Ct,42.50.Dv,03.67.Mn 



I. INTRODUCTION 

EfRcient entanglement generation, distribution and 
storage are the bàsic requirements for development of 
successful quantum information technologies, including 
proposals for realizing quantum information networks 
They comprise entanglement generation by Raman 
scattering, storage in collective polarization of atòmic 
ensembles and release by anti-Stokes Raman scattering. 
Such a network can be used, among other things, for 
a quantum cryptography Q and quantum teleportation 
Those protocols deal with spin or polarization en- 
tanglement. Another important possibility is the entan- 
glement of continuous variables, such as light quadrature 
amplitudes or molecular vibrational coordinates. Stokes 
Raman scattering [5, 6] generates entanglement between 
light and àtoms, and this in turn can be used for entan- 
gling distant atòmic ensembles @, El El- Such an entan- 
glement can be exploited for continuous-variable quan- 
tum cryptography ° r teleportation \x&- In the Stokes 
or anti-Stokes configuration, Raman scattering can be 
used for transferringa quantum state of light to àtoms 
or vice versa (l2L IÏ3I. IÏ4L l·L5l | . Already early works demon- 
strated phase memory of atòmic polarization induced by 
Raman scattering and its robustness against decoherence 
[ïfjj . thus a Raman medium offers the possibility to real- 
ize long-lived continuous-variable quantum memory. 

Parallel to Raman scattering techniques, the phe- 
nomenon of Electromagnetically Induced Transparency 
(EIT) is being extensively studied. In this completely 



resonant situation, weak probe light can be slowed down 
and stopped in a form of localized atòmic polarization by 
controlling a strongpump [ÏtI flsl Iïfl| . In this form of a 
quantum memory [2fJ|, the light bandwidth is narrower, 
compared to the off-resonant Raman interaction, due to 
the resonant character of this phenomenon Q] . 

In this paper we focus on the precise quantum de- 
scription of the generation of entanglement between light 
and an atòmic ensemble by spontaneous and stimulated 
Stokes Raman scattering, followed by the readout of the 
atòmic polarization state onto an anti-Stokes field, as de- 
picted in FigQ]and[21 Following Raymer Q, we will use 
an input-output formalism to describe these interactions. 
As was already conjectured there and proved approxi- 
mately by numerical means, the Raman scattering pro- 
cess can be decomposed into a discrete set of two-mode 
entanglement processes. Each process is described by 
a Bogoliubov two-mode squeezing transformation, which 
linearly mixes photon annihilation and creation opera- 
tors. Each elementary process squeezes the phase sum 
of excitation amplitudes of optical and atomic-ensemble 
collective bosonic modes. Below we shall present a proof 
of this decomposition, an extension of the Bloch-Messiah 
reduction 21], and demonstrate that it is applicable as 
long as the input-output relations are linear relations, 
involving mixing creation and annihilation operators. 

We also consider the readout of the stored atòmic po- 
larization using the anti-Stokes scattering process. We 
prové that the readout process can also be decomposed 
into múltiple mode pairs, each pair undergoing indepen- 
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dent evolution analogous to a beam-splitter transforma- 
tion. For this, we derive a canonical form of an all-linear 
unitary transforation acting on a bipartite system. 

Throughout this paper we keep in mind the idea of an 
experiment with Raman scattering in a thermal vapor of 
87 Rb, pumped by highly detuned, sub-nanosecond làser 
pulses. Our idealization comprises neglecting the influ- 
ence of sidewards spontaneous emission, as well as atòmic 
polarization damping, atom loss and density fluctua- 
tions. We will also restrict ourselves to a one-dimensional 
model, which is justified approximately in the case when 
the Fresnel number of the pump beam equals unity or 
less jfj]. However, we will include the effects of dispersion 
in the vapor, which introduces appreciable group-velocity 
difference between the pump and scattered wave under 
realistic experimental conditions. We find that disper- 
sion signiflcantly increases the number of independent 
modes excited in the Raman scattering process, since 
the Stokes photons born in different parts of the medium 
become distinguishable by their delay respective to the 
pump pulse. Dispersion also can reduce the emciency of 
the readout process. 

Our work stems from early studies of stimulated Ra- 
man scattering 0, 0, an d we adopt the formalism 
of these works. We will also draw from early results 
on quantum energy [23L l2~3 | and phase j2^| statistics of 
Stokes pulses. 

In Sec.[ïï]we give theoretical foundations of the decom- 
position used to analyze the structure of squeezing gener- 
ated in the Stokes scattering process. The specific case of 
pulsed pump fields and nonzero dispersion, which needs 
to be solved numerically, is addressed in Sec. IIIII The 
consequences of the multimode character of the output 
field for photon-count statistics are studied in Sec. IIVI 
The problem of reading out of the atòmic polarization 
by anti-Stokes scattering is discussed in Sec. 
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FIG. 2: Pictorial representation of the Stokes and subsequent 
anti-Stokes scattering processes (see Sec. in the labora- 
tory reference frame. The vertical stripe represents the celi in 
which the nonlinear medium is contained, the shaded paral- 
lelograms are the spatio-temporal interaction boundaries, the 
grey (green in color) stripe represents Stokes pump field enve- 
lope A p (z,t), while the dark grey (red) stripe represents anti- 
Stokes pump field envelope A' p (z,t). Letters with hats are 
the annihilation operators of the input and output quantum 
fields, defined at the respective boundaries of the interaction 
regions. 



II. PAIR-WISE ENTANGLEMENT OPERATION 
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FIG. 1: Relevant atòmic levels in the Stokes and anti-Stokes 
(right panel) Raman scattering. Initially àtoms reside in state 
1). During the first scattering process (left panel) they are 
driven by strong pump field at a frequency u> p , which creates 
Stokes optical field à at a frequency lü 3 and atòmic polariza- 
tion b. Subsequently, the àtoms are subjected to a strong field 
at a frequency uj p i which erases atòmic polarization we will 
denote by d shifting the excitation into anti-Stokes optical 
field c at a frequency ui a (right panel). 



We begin developing our model by idealizing the pro- 
cess of stimulated Raman scattering as a three-wave mix- 
ing and restricting ourselves to one dimension This 
is justified provided that atòmic saturation and pump 
field depletion can be neglected while the pump Rayleigh 
range is comparable to or shorter than the medium length 

@ 

We formulate the equations in the moving reference 
frame of the Stokes wave 

z-L/2 

t = ÍLAB — (1) 

^gr,s 

where ílab is the time in the laboratory frame, v gItS 
is the Stokes group velocity, while L is the medium 
length. According to the above expression, the labo- 
ratory and Stokes reference frames meet in the middle 



of the medium. Atòmic polarization (or spin is de- 
scribed by an operator b(z,t) which removes one atom's 
worth of internal excitation from the atòmic ensemble 
localized in a thin-slice volume around a point in space- 
time in the medium [jjj]: 



b(z,t) = i?— V |l a )(3 a | 

71 ' 



1/2 



J(ui p -Lj s )t L·AB -i(k p -k s )z 



(2) 



where p is linear atom density (àtoms per mm), n is the 
number of àtoms in a slice, indexed by a, |1) and |3) de- 
note the Raman-coupled atòmic levels, ui p and uj s are the 
pump and Stokes angular freqüències, while k p and k s are 
their wavevectors (see Fig.QJ. As long as the fraction of 
excited àtoms is small, b(z,t) is approximately a bosonic 
operator, [b(z, í), 6(z', t)>] sa 8{z — z'). The Stokes field 
is described by an envelope annihilation operator at each 
point in space along the pump beam a(z,t): 



à( Z ,t) = / düJ e -Í("-".)*L 



z à(z, Lü) 



(3) 



It is a bosonic operator, [à(z, t), à(z, í')t] = S(t — t'). Fi- 
nally, the pump field is given by its envelope A p (z, t) in 
the reference frame of the Stokes pulse, scaled to unity, 
defined by 



E p (t - {z - L/2)A0) 



E, 



= A p (z, t)e 



—iüjpt-\-ik p z 



(4) 



p,max 



where E p {t) is the pump elèctric field while A/3 is the 
difference of the inverse group velocities of the pump and 
Stokes pulses (measured in picoseconds per milimeter) 



A(3 = 



1 



(5) 



where v gTtP is the pump group velocity. If the Stokes field 
is faster than the pump, then A/3 is negative. 

Using these definitions, the Raman scattering process 
is cast into the following pair of coupled equations 



dà(z, t) 
dz 



= 9oA p (z,t)b\z,t) 



= g A p (z,t)iï(z,t) 



(6a) 
(6b) 



where go is the coupling constant (measured in 
(mmps) -1 / 2 ) defined as: 



90 = P 1 ' 2 



flüJ s 

2e c 



(7) 



where 5" is the pump beam cross-section while k,\ is a 
coupling constant defined in Q. Since the pump and 
Stokes pulses are narrowband, higher order dispersion 
effects, possibly distorting the pulses, can be neglected. 
On the other hand, in the atòmic vapor one finds that 
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FIG. 3: Pictorial representation of the Stokes scattering pro- 
cess in the Stokes reference frame, defined in Eq. The 
medium entangles light field à with àtoms b. The local entan- 
glement amplitude increases with a (fixed) local pump field 
strength A p (z,t), represented in the picture by a grey (green 
in color) stripe. Horizontal stripe corresponds to equal group 
velocities of pump and Stokes, A/3 = 0. 



the frequency difference between the interacting waves is 
large enough for significant group velocity differences to 
be observed. A typical subnanosecond pump pulse can 
be delayed by several times its duration with respect to 
the Stokes pulse, thus the presence of nonzero A/3 cannot 
be neglected [26l |, 

We proceed with solving Eq. 10 by integrating over 
a finite spatial-temporal region in which the interaction 
takes place, thus formulating the problem as a quan- 
tum scattering problem. This approach, pictorially rep- 
resented in Fig.|21requires introduction of a suitable scat- 
tering matrix, which in the case of Stokes scattering is 
composed of the Green functions of Eq. ©. We choose 
to introduce them over a space-time region < z < L 
and —T < t < T, in the following way: 



à(L,t) = / dt'C a (t,t')à(0,t') 

J-T 

+ í dz'S a (t,z'p{z',-T) 



b(z,T)= dz'C b (z,z')b(z',-T) 



+ / dt'S b (z,t')à^(0,t') 



(8b) 



We can simplify the above input-output relations us- 
ing an extended Bloch-Messiah reduction theorem [2lj | . 
which we derive in the AppendixEI This theorem states 
that since the Green functions for this problem must 
describe a unitary operator transformation, they pos- 
sess common singular vectors and related singular vàlues 
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cosh ( n and sinh : 






C a {t,t') =EcoshC„ 

II 


•^*(í)^ in) (í') 


(9a) 


C b (z,z) = E cosn Cn 

n 


■ <f>í out) *(z)4: n) ^) 


(9b) 


S a (t,z') = E sinh Ç n 

n 


^ n out) *(t)^ n >(z') 


(9c) 


S b (z,t') = E sinh Ç ra 




(9d) 



where ( n are real numbers while ipn (t), 4>n Ut \z), 

ipn^ (í) and Çn (z) are complex- valued functions which 
form separate orthonormal bases for the output and in- 
put, light and atòmic mode-function spaces respectively 



dt^ n >(t)-^\t) = 5 n , 



dzcj>t ) *{z)-4>t ) {z) = 5 n , 

r 

dz<f>i out >(z)-cf>(™V(z)=6 n , 



(10) 



We insert the Green functions in the form given in Eq. @ 
into Eq. |JHJ. Then we express input and output field an- 
nihilation operators à(0, t), b(z, — T) by mode expansions 
and define output operators à(L,t) and b(z, T) analo- 
gously 



o(0,í) 


n 


v4 in) *(*) 


à(L,t) 


n 


)^ oui) *(í) 


b{z,-T) 


n 




b(z,T) 


= E^ 4 


fà ut >(z), 



(11) 



where we have introduced operators àn , àn Ut \ bn^ 
and. bn . Thanks to orthonormality of the respective 
mode functions expressed in Eq. ifÏQl) . these operators 
obey canonical commutational relations: 

[àí out) , aÍT í)f ] - 5 nm [bir t} , bÍT m ] = Snm (12) 
and the operators with superscript (in) obey analogous 
commutational relations. Therefore 
and bn are bosonic operators, which annihilate excita- 
tions in respective light and atòmic modes. Using Eq. ©, 
(1 1 01) and lli lli . Eq. © can be cast into the following form: 

ü {out) = à^coshCn + òl m)t sinhC„ (13a) 
y(out) = ^coshCn + ò^^sinhC,,. (13b) 



This expresses independent two-mode squeezing pro- 
cesses, and is the main result of this part of our study. 
Each of these elementary processes squeezes the opti- 
cal field in a particular temporal mode (í) with the 
atòmic polarization in a particular longitudinal spatial 
mode 4>n 7 ^ (z), producing fields in optical temporal modes 
^n° ut) (t) that are pair-wise entangled with the internal 
states of àtoms described by collective spatial modes 

4>n Wt \z). The squeezing parameter equals Cn ■ The noise 
in the sum and difference of quadratures scales by a factor 
of exp(£„), corresponding to the mean number of excita- 



tions (n n ) = (à 



(cmí)t - (out) 



> = (&: 



On 



sinh Cn- 



Note that for a pump pulse of a given shape, charac- 
terized by duration t p , for instance Gaussian: 



A p (z, t) = exp 



-2 ln 2 



t- (z-L/2)A0 



, (14) 



the time and space in equations © can be reduced 
to a dimensionless coordinates by measuring time in 
units of pump pulse duration t p and space in units 
of distance over which temporal Stokes walk-off equals 
pump duration T p /A/3. This way the pump envelope 
A p (A(3z/Tp,t/T p ) becomes a fixed function. Note that 
for T >• t p the Green functions do not depend on T 
since there is no longer any interaction taking place for 
large t. Thus the Green functions defined in Eq. © can 
depend only on medium length, coupling strength, and 
pulse shape. Both can be brought into a form of two 
dimensionless parameters: 



LA/3 



(15a) 
(15b) 



The first is the temporal walkoff between pump and 
Stokes field, measured in pump pulse durations, while 
the second measures the squeezing strength and is pro- 
portional to the squeezing parameter in the case of zero 
dispersion 0, 0|. Therefore there is a congruence be- 
tween the solutions and mode functions for the Stokes 
scattering problem, as long as the pump pulse shapes are 
congruent and the vàlues of the parameters T and A are 
the same. 

Finally, let us note that the characteristic output 
modes ipn° (t) and </4 ouí · ) (t) we obtained with the help of 
the Bloch-Messiah reduction are identical to previo uslv 
obtained atòmic and field coherence modes 0, 0, l27| . 
The latter were introduced as the eigenmodes of the first- 
order coherence functions for àtoms (b(z, T)w(z', T)) and 
field (à(L, t)à'(L, t')). Calculating these functions explic- 
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itly using equations 10 we obtain: 

(à(L,t)a\L,t')) = f L dz'S a (t,z')S:(t',z') 
Jo 

= J2(nn)€ OUt) *(t)€° Ut) (t') (16) 

n 

(b(z,Tp{z',T))= í dt'S b (z,t')SZ(z',t') 

J-T 

=Y,^)4r t) *{t)<l· [ r t) {t i ) (i?) 

n 

as in 0, 0| . In also the characteristic input modes 
were introduced as those evolving into a particular pair 
of output modes, which is another way of obtaining the 
decomposition given in Eq. ©. 



III. NUMERICAL RESULTS 
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FIG. 4: Mean number of Stokes photons as a function of cou- 
pling strength go [(mm ps)~ 1,/2 ]. Three lines are the results of 
computations for different A/3 equal (sòlid line), — lOps/mm 
(dashed line) and —30 ps/mm (dotted line). 



To illustrate the meaning of the above results and pro- 
vide the decomposition into independent entanglement 
processes in some specific and realistic cases, we have 
numerically found Green functions deflned in Eq. ©. 
Calculations were carried out for a medium of length 
L = 75 mm for various vàlues of group velocity differ- 
ence A/3, coupling strength go and Gaussian pump pulse 
of FWHM duration t p = 200 ps as deflned in Eq. ifLlj) . 

As the equations of motion 10 are linear, they are iden- 
tical to the classical equations for evolution of the Stokes 
fleld and the atòmic polarization in the case of stimulated 
Raman scattering, in which annihilation operators à(z, t) 
and b(z,t) are replaced by classical fleld a(z,t) and po- 
larization amplitudes f3(z,t). Therefore we can compute 
matrix approximations to the Green functions C (í, i'), 
Cb{z,z'), S a (t,z') and Sb(z,t') using Standard methods 
developed in nonlinear òptics. To accomplish this, equa- 
tions © were solved numerically for a complete orthonor- 
mal set of initial conditions. First we put a(0,t) = for 
every t except for a single point on a computational grid 
t = t' where a(0,t') = 1, and f3(z,—T) — for every z. 
Performing finite integration steps over the space-time 
grid of 200x200 points where the interaction occurs, we 
computed a(L,t) and (3(z,T). For the initial conditions 
chosen, they are equal to C a (t,t') and Sb(z,t'), respec- 
tively. Repeating this calculation for each t', we con- 
structed a matrix approximation of these Green func- 
tions. Then analogously starting from a(0, t) — 0, and 
j3(z, — T) = except for a single point on a computational 
grid where @(z', —T) = 1, we have constructed Cb(z, z'), 
S a (t,z'). Next we numerically performed the singular 
value decomposition of these matrices to obtain the 
singular vàlues sinh£„ and cosh£„ as well as associated 
vectors ^(t), ^(t) and ^ n ° ut) (t) as de- 

flned in Eq. ©. We have verifled that our simulation is 
accurate by repeating it on a grid two times flner in each 
dimension. The differences between the results were of 



the order of precision with which singular value decom- 
position can be computed. We have also verifled that our 
numerical calculation reconstructs the analytical results 
for square pump pulse without dispersion 9j. 

Numerical calculations confirm that the total number 
of photons 

(Ntot) = X>»> = sinh2 ( 18 ) 

n n 

is, to a good approximation, an exponential function of 
coupling go, within the range plotted in Fig.QJ However, 
we flnd a strong dependence of (iV tot ) on the group veloc- 
ity difference between pump and Stokes pulses A/3. With 
increasing A/3 the total photon number (N tot ) quickly de- 
creases. Using results displayed in Fig.Q]we have found, 
for each A/3, the value of coupling go for which the total 
number of photons equals (N tot ) = 10 6 . These vàlues of 
coupling go will be used in the following examples. 

In Fig. we plot the mean number of photons in each 
mode (n n ) = sinh 2 £„. In the case of no dispersion 
A/3 = 0, the occupation of the first mode dominates 
by a factor of 10 3 over the next, as known previously 
[23,|2J|. Strong domination, however, quickly comes to 
an end with increasing dispersion. For A/3 = —30 ps/mm 
more than a dozen modes have occupancy equal within 
an order of magnitude. 

The physical reason behind this result can be traced in 
Fig. El where we depict a space-time region where the in- 
teraction occurs. We mark the pumped area as the dark 
grey (green) region. In the case of zero group velocity 
difference, A/3 = 0, the first few Stokes photons propa- 
gate through the pumped region, stimulating subsequent 
scattering acts and inducing coherence between atòmic 
polarization in various regions along the celi. The pro- 
cess is almost single mode On the other hand, 

in the case of high group velocity difference, the spatial- 
temporal interaction region can be divided into several 
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FIG. 5: Mean number of photons in each of characteristic out- 
put modes (n n ) for 15 most-occupied modes. The histograms 
has been made for coupling yielding average total number of 
Stokes photons (N to t) = 10 6 for different A/3 equal (black 
bar), —10 ps/mm (grey bars) and —30 ps/mm (empty bars). 



is time delayed relative to pump peak [23. The modes 
resemble in shape the eigenmodes of a harmònic oscilla- 
tor. They are also strongly afïected by the group velocity 
difference. With increasing A/3, the duration of the fun- 
damental mode increases while the time delay relative to 
the pump pulse in the middle of the medium decreases. 
Also for large A/3 all the characteristic modes have sig- 
nificant intensity over the same period of time, spanning 
approximately one pump duration plus temporal walkoff 
between pump and the Stokes r p + LA/3 = t p (1 + A). 

Note that in the limit of zero dispersion, A/3 = 0, the 
Green functions calculated numerically can be found in a 
closed form using Laplace-transform techniques Q, even 
for the full three-dimensional case @. The results pre- 
sented here for zero dispersion are recovered by applying 
numerical diagonalization directly to those functions (§]. 
Finally let us note that the output characteristic modes 
can be used to reconstruct the statistics of fluctuating 
Stokes pulse shapes (2^ |. 




FIG. 6: Pictorial representation of the Stokes scattering for 
zero group velocity difference A/3 (left) and large negative 
A/3 -C — Tp/L (right). The grey box represents the spatio- 
temporal region of interaction, the dark grey (green in color) 
stripe is the pumped region, the horizontal (red) arrows are 
the Stokes emission, while the vertical (blue) arrows represent 
atòmic polarization. The picture has been drawn in the Stokes 
reference frame. 



regions, the number of such regions being roughly the 
ratio of the total temporal walkoff to the pump duration 
A = LA/3/t p . These regions are weakly coupled, since 
the photons scattered in the distinct parts of the celi do 
not have a chance to meet each other. Thus several un- 
correlated temporal modes take part in the process and 
become significantly occupied. 

This latter phenomenon complicates the possible appli- 
cations. For nonzero dispersion one must afford means, 
such as for example balanced homodyne detection, for 
separating distinct output modes in order to extract a 
pure squeezed state of light and matter. Whether in the 
single or multimode squeezing regime, an efficient detec- 
tion of Stokes light is necessary for quantum communi- 
cation and computation. This is possible by means of 
balanced homodyne detection using a local oscillator in 
the shape of the characteristic Stokes mode. In Fig.[3we 
plot the temporal shape of the most occupied output field 
modes for the Gaussian pump centered at t = defined 
in . The dominant mode has a single peak, which 



IV. PHOTON COUNT STATISTICS 

Using the results of numerical calculations, we were 
able to compute the photon-count statistics for the Stokes 
pulse. Let us recali that one mode of the pair compris- 
ing a two-mode squeezed state exhibits thermal statis- 
tics. This is the case of any particular output field mode 
ipn (t) by itself. Also, each of the output field modes is 
independent, since they evolve separately from vacuum. 
Thus the Stokes pulse has multimode thermal statistics 
El: 



P(n) = i 



-, — ihr> 



dk 



n„(i-<*<a»>) 



(19) 



We plot p(n) in Fig.|Slfor total photon number (N tot ) = 
10 6 and various vàlues of group velocity difference A/3. 

Following the analogy with thermal light, we have also 
computed and plot in Fig. the equivalent number of 
modes Ai. It is defined as a number of independent 
thermal modes of equal intensity that measured together 
reproduce the Stokes fluctuation-to-mean ratio: 



M 



(N t , 



Er, 



(20) 



As can be seen in Fig. El the equivalent mode number Ai 
depends mostly on A/3, which is almost independent of 
coupling strength g . This is consistent with the intuitive 
way of understanding Stokes scattering we propose in 
Fig. El 



V. ATÒMIC POLARIZATION READOUT 

As shown in Section[ïï] the optical Stokes pulse con- 
tains only one part of the squeezed state. Alone it has 
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FIG. 8: Photon count probability distribution p(n) of Stokes 
light pulses. Three lines are the results of computations for 
different A/3 equal (sòlid line), —10 ps/mm (dashed line) 
and —30 ps/mm (dotted line). 
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FIG. 9: Equivalent number of independent thermal modes Ai 
as a function of coupling strength go [(mm ps) -1 / 2 ]. Three 
lines are the results of computations for different A/3 equal 
(sòlid line), —10 ps/mm (dashed line) and —30 ps/mm (dotted 
line). 



to Eqs. © in the anti-Stokes reference frame 



thermal statistics and possesses no distinct quantum fea- 
tures. Entanglement can be observed only if the field 
counterpart — the atòmic polarization — is acted upon. 
A convenient to way process the atòmic b(z) field is to 
translate its statistics into the optical field. This is done 
in the anti-Stokes scattering process, during which a por- 
tion of atòmic polarization gives rise to an anti-Stokes 
pulse. The equations governing this process are similar 



= g' A' p (z,t)d(z,t) (21a) 

= - g '*A';(z,t)c(z,t) (21b) 

where c(z,t) is the anti-Stokes field annihilation opera- 
tor, defined by an equation analogous to (pj, d(z,t) is 
the atòmic polarization operator, defined by an equation 
analogous to (J2J , g' is the readout coupling defined anal- 
ogously to |0, while A' p (z,t) is the readout pump field 
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envelope. The definition of the A' [z, t) is analogous to 
that of the Stokes pump field envelope A p (z,t) given in 
Eq. but with a different peak amplitude and a dif- 
ferent group velocity difference, which we will designate 
by A/3': 



A/3' 



u gr,a 



u gr,P 



(22) 



where v gI . ja is the group velocity of the anti-Stokes wave, 
while v gl p is the group velocity of the readout pump. 

Again, we formulate the readout problem as a quan- 
tum scattering problem, analogously to the situation il- 
lustrated in Fig. This time the scattering matrix is 
composed of four Green functions for the anti-Stokes 
scattering equations 12111 . We define them analogously 
to those deflned in Eq. <|HJl for the Stokes scattering: 



c(L,t) 



d[z,T) 



dt'C c (t,t')c(0,t') 



+ / dz'S c (t,z')d(z',-T) (23a) 



(23b) 



dz'C d {z,z')d{z',-T) 



— T 



dt'S d (z,t r )c(0,t') 



Here c(0, t) is the optical field operator (usually repre- 
senting the vacuum) entering the vapor during readout, 
and d(z, — T) is the atòmic ensemble operator just prior 
to the readout. 

We can cast these relations into múltiple, independent 
beamsplitter transformations (see Appendix [Bj . Since 
the Green functions for this problem must describe a uni- 
tary transformation, they possess common singular vec- 
tors and related singular vàlues, which we parameterize 

by T] n - 



C c (t,t') 


n 


-vn-Mr t] *(t) *l in) (*') 


(24a) 


Cd(z, z) 


=EVi- 

n 


- Vn -^ n ° ut >(z)^: n \z') 


(24b) 


S c (t, z') 






(24c) 


S d (z,t') 




■&° ut) *{z)^ n) {t'). 


(24d) 



The readout modes ¥ n out \t), <S> ( n Ut \z), *£ n) (t) and 
&n n \z) form orthonormal bases in the intput/output, 



light and atòmic modefunction spaces: 

f-T 



dt^ n >(t)-¥^ n \t)^5„ 



T 
L 



-T 
L 



(25) 



dz&: n >(z)-$W(z) = ó n , m 
dt^° ut >(t)-&™ t \t) = ó nim 

dz^ n out >(z)-^ t \z)=6^ m . 

We use these functions as mode bases in respective spaces 
and decompose the annihilation operators c(0, í), c(L,í), 
d(z,~T) and d(z,T): 

g(0,í) = (26a) 

n 

c(L,t) = Y / à { r t) ^r t) *(t) (26b) 

n 

d(z^T) = Y / d i : n) ^: n) *(z) (26c) 

n 

d(z, T) = ^2 d ( n ut) <S>í° ut) *(z). (26d) 



and dn^ 

~Í0Ut) j(out 

tors c„ and d n 



Where we have introduced annihilation operators of the 
input and output characteristic modes Cn , <%n , Cn 

Again, it turns out that the output opera- 
are independent bosonic operators. 
Same is true for input operators cl™'* and dn"^ . 

These operators annihilate excitations of modes that 
undergo transformations analogous to a beamsplitter 
transformation: 



d ( n out) = V^üà^ - 



s("0 



'In L-n 



(27a) 
(27b) 



The real parameters rj n have the meaning of the readout 
efficiency. Modes with r\ n = 1 are completely translated 
from the atòmic field into the optical field while spatial 
modes with r\ n = remain unaffected by the anti-Stokes 
scattering process. 

We have numerically found forms for the Green func- 
tions defined in Eq. (1 2 - > I for similar conditions used to 
calculate the Raman process in earlier sections. Results 
are first presented for the case of a medium of length 
L = 75 mm with zero group velocity difference A/3', 
various coupling strengths g' Q and Gaussian pump pulse 
A! p {z,t) of FWHM duration t' p = 200 ps given in Eo. ltÏ4|l . 
In Fig.^lwe plot the calculated r\ n . It can be seen that 
above a certain value of coupling g' , many modes have 
unit readout efficiency 7?„ = 1, i.e. are completely trans- 
lated into the optical field. Those modes are degenerate, 
thus within a subspace spanned by their mode functions 
we can change the basis, i.e. the modal functions, with- 
out affecting relations l|27JI . 
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FIG. 10: Readout efficiency r\ n for the most efficiently read 
modes in anti-Stokes scattering process for A/3' = and 
coupling g' equal 0.01, 0.02, 0.05, 0.1 (psmm) -1 ' 2 , for black, 
grey and empty bars respectively. 



Let us focus on the case that the readout takes place 
just following the Raman scattering process described in 
earlier sections; in this case d(z, —T) is given by the out- 
put operator b(z,T) in Eq. JHJ|: 



d(z, -T) 



b(z,T) 



(28) 



where Afc is the wavevector mismatch in the time-delayed 
four-wave mixing process: 



Afc 



(29) 



where subscript p, p', s, a denote Stokes pump and anti- 
Stokes pump, Stokes and anti-Stokes signal waves respec- 
tively. For efHcient readout the phase-matching condition 
Afc = must be satisfied. Below we will assume perfect 
phasematching. It is possible to satisfy this by proper 
choice of pump freqüències. 

Using the expansions 11 lli and i'K'ú) we find that the 

characteristic anti-Stokes input operators dn^ are re- 
lated to the characteristic Stokes output operators bn Vt ^ 
by an unitary transformation expressing the change of 
basis from <fin (z) to (z): 



j(in) _ V^tt Uout) 
"m / , u mn"n 



where the transformation coeficients are: 



dz4r th (z)e- 



■iAkz 



■ 4> 



(in) 



(z). 



(30) 



(31) 



Let us consider the readout of the flrst characteristic 



atòmic mode of the Stokes scattering 



,(out) 



(z) whose 



quantum statistics, before the readout, is contained in 
the operator ^ ouí ) . It can be decomposed in the charac- 
teristic anti-Stokes input mode base $n (z), according 



to the above formulas. The readout has unit efficiency if 
(z) decomposes only into readout 



the atòmic mode <^°"^ 



characteristic input modes Qn n \z) for which i) m = 1. In 
other words, if for all m for which U m i ^ the read- 
out efficiencies rj m equal one, then the readout has unit 
efficiency. 

There is however a more straightforward way to ana- 
lyze a readout of an excitation in a particular mode. As- 
suming perfect phase matching, Afc = 0, so d(z, — T) = 
b(z, T) we substitute input atòmic operator d(z, — T) in 
Eq. (l23l) by b(z, T) in an expanded form given by Eq. itïïll : 

c(L, t)=Y, <(0n Ut) + [ dt'C c (t, t')c(0, t') (32a) 
d(z,T) = J2e* n (z)b^ - i dt'S d (z,t')c(0,t') (32b) 

n ~T 

where we have introduced mode functions 

a* n (t)= [ L dz'S c (t,zi)^*(z') (33a) 
Jo 

s* n (z)=í dz'C d {z,z')4> { ° ut) *^')- (33b) 



The equations l(32|l give yet another modal decomposition 
of the output atòmic and light fields. However this time 
the modal functions: cr*(í), C c (í, í'), e„(z) and Sd{z,t') 
form othonormal bases in the joint atom-light space: 

dt(j* n (t)cF m (t) + / dze* n (z)e m (z) = S mn (34a) 
-t Jo 

T r-L 

dta* n (t)C c (t,t')- dze; i (z)S d (z,t r ) = 0. (34b) 



This is a consequence of unitarity of the relations ((22 
(see also Annendix lB)l . 

Usually we want statistics of some particular atòmic 
mode, say </>i°"^(z), to be completely transferred into 
an optical field mode. This requirement is equivalent 
to requesting e\{z) — for every z, so that d(z,T) in 
Eq. l|32bjl contains only zero-point, vaccum noise repre- 
sented by c(0,t'). The output light mode eri(t) can be 
computed in advance and detection apparatus tuned to 
it. Note, that if the transfer is complete the z integrals 
containing S\{z) vanish in Eq. 1341 . It follows that all 
other output field modes must be orthogonal to oi(t). In 
other words no crosstalk occurs, i.e. the polarization in 
any other atòmic mode does not give rise to any photons 
in the optical field mode <J\{t) we are presumably tuned 
to. This does not prevent emission into other optical 
modes, however. 

Therefore the readout problem can be reduced to a 
question of how to prepare the readout pump pulse in 
such a way that the atòmic polarization of interest, say in 
mode 4>^ ut \z), is completely transferred into the optical 
field. 
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FIG. 11: Praction of residual atòmic polarization 1 — rj not 
transferred to an optical field after the anti-Stokes scattering 
process as a function of the readout coupling go [(ps mm) -1 / 2 ]. 
Results for three group velocity differences A/3' equal (sòlid 
line), —10 ps/mm (dashed line) and —30 ps/mm (dotted line) 
are displayed. 



We investigated this by simulating readout of an 



atòmic mode 



, (out) 



(z) in an anti-Stokes interaction. Cal- 



culations were carried out for a medium of length L = 
75 mm for various vàlues of group velocity differences, 
assuming A/3 = A/3'. For each particular value of this 
parameter, we have first computed the mode function 
of the fundamental atòmic output mode <j)f ut ^ (z) in a 
way described in Sec. IIIII with the Stokes coupling such 
that the total mean number of excitations (N tot ) = 10 6 . 
Next we have simulated the anti-Stokes scattering of this 
particular atòmic mode for various anti-Stokes coupling 
strengths g' Q and Gaussian pump A' p (z,t) of FWHM du- 
ration r p = 200 ps given in Eq. lUljl . We have quantified 
the quality of the readout by computing the residual frac- 
tion of atòmic polarization after the readout 



1 



n : 



dz\ei(z)Y 



(35) 



The results are shown in Fig. El O ne can see that for 
zero group velocity difference 1 — rj rapidly falls with in- 
creasing coupling g' Q as predicted in The readout is 
virtually perfect for g' > 0.04 (psmm)^ 1 / 2 ). However 
for nonzero A/3' the residual atòmic polarization is of 
the order of at least a percent. The readout inefEciency 
l—i] becomes a periòdic function of coupling g' . The 
period increases with increasing A/3', and the particular 
optimal vàlues of g' depend on the parameters of the 
readout process as well as the initial atòmic excitation 
mode function (j)^ ut \z). 

In Fig.^Jwe plot the output field modes <J\(t) obtained 
in the readout simulations. In case of homodyne detec- 
tion of the output of the anti-Stokes readout process, 
those would be the optimal local oscillator shapes. For 



zero group velocity difference A/3' = we observe that 
the output field mode becomes shorter with increasing 
g' . For nonzero A/3' the mode shapes 01 (i) correspond- 
ing to the subsequent minima of the function 1 — r](g' ) 
depicted in the Fig.EJshow increasing number of nodes. 

The physics behind the readout process can be in- 
tuitively understood if one notes that integrating the 
propagation equations 12H over a spatiotemporal region 
2Az x 2Aí where the fields are almost constant gives a 
simple beamsplitter-like relation between input and out- 
put fields: 



c(z + Az, t) =c(z - Az, t)y/ï - I7I 2 + jd(z, t - Ai) 



d(z, t + At) =d{z, t - Aí)VÏ - bf - l*c{z - Az, t) 



At 
Ai 



where 7 = 2g' A' p (z, í)V AzAí. Therefore we can visual- 
ize the scattering process as a square network of virtual 
beamsplitters, each having different reflectivity ^7. In 
Fig-E3 we represent the pumping intensity A' p (z,t) and 
the beams representing atòmic polarization and optical 
field, which propagate in above mentioned network. 

For zero group velocity difference and a nearly expo- 
nential input atòmic mode </)^ ont ' ) (z) (see Fig.Gk) the anti- 
Stokes field grows with z sweeping the major part of the 
atòmic polarization out of the celi. At each z along the 
pump beam, where the virtual beamsplitters mix atòmic 
polarization and anti-Stokes light the condition of sup- 
pressing the atòmic output d(z, í)\/l — I7I 2 = 7*c(z, t) is 
met, since the anti-Stokes and atòmic polarization both 
grow exponentially with z. For small group velocity dif- 
ference A/3' = —10 ps/mm, the input atòmic polarization 
has a maximum inside the celi and the anti-Stokes scat- 
tering process first shifts the atòmic excitation into the 
field, but at some point the process is partially reversed 
and the atòmic polarization, albeit shifted towards the 
end of the celi, partially remains within it. The atòmic 
input mode shape is not proper to achieve suppression 
of the atòmic output polarization which occurred in the 
case of the of vanishing dispersion. Finally, for the largest 
group velocity difference A/3' = —30 ps/mm, the initial 
atòmic mode shape plays a minor role, since the scat- 
tering processes in various parts of the celi are nearly 
independent from each other, because there is no light 
connecting them. This is why the residual fraction of the 
atòmic polarization can be smaller for larger A/3' as found 
in the numerical calculations and depicted in Fig. El 

We explain the oscillations of readout efficiency 1 — 77 as 
a function of coupling g' found in numerical calculations 
and plotted in Fig. El m the following way. We investi- 
gate in detail the scattering in a small part of the pumped 
region, shown in Fig. I13H . As soon as the atòmic polar- 
ization reaches the pumped region, it is converted into an 
anti-Stokes wave within a small period of time. Released 
light propagates along the horizontal arrow, however soon 
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FIG. 12: Readout anti-Stokes light modes <ri(í) for fundamental atòmic Stokes modes 0i (z) for various group velocity 
diíïerence A/3 = A/3' and coupling strengths corresponding to readout efficiency peaks for nonzero A/3' (compare Fig. I1H . 
Sòlid, dashed and dotted lines correspond respectively to subsequent vàlues of g' , listed below each panel. 



it is converted back into atòmic polarization. The exci- 
tation exchange continues until the atòmic polarization 
or light leaves the pumped region. Since the period of 
oscillations changes with the pumping intensity, we ob- 
serve in the end either excitation form alternately. In 
particular for some pumping intensities the oscillations 
stop just when the excitations are stored in light, which 
corresponds to the maximum readout efficiency 77. This 
scenario is confirmed by the numerical simulations for 
high pumping intensity. It can be also obtained from an 
analytical solution of Eq. 112111 rewritten in the pump ref- 
erence frame under assumption that the quantum fields 




(c) (d) 

FIG. 13: Pictorial representation of the anti-Stokes scattering 
for zero group velocity difference A/3' (a), medium value of 
A/3' (b) and high negative A/3' <g — r p /L (c). The grey box 
represents the spatio-temporal region of interaction, the dark 
grey (red in color) stripe is the pumped region, the horizontal 
(green) arrows are the anti-Stokes emission, while the vertical 
(blue) arrows represent atòmic polarization. The picture is 
drawn in the anti-Stokes reference frame. Panel (d) gives 
insight into oscillatory exchange of excitation between light 
and àtoms. 



do not depend on z in this reference frame. 



VI. CONCLUSION 

In summary, we have analyzed the Stokes and anti- 
Stokes scattering processes under transient pumping con- 
ditions and in presence of group velocity difference be- 
tween the interacting waves in the non-saturated regime. 
Using Bloch-Messiah reduction |2ÏJ we confirm previous 
conjecture 0| that the Stokes scattering process can be 
decomposed into múltiple independent squeezers. Each 
of them entangles àtoms and field, both in characteristic 
input modes ipn (í) and 4>n (z), into a pair of quantum- 
correlated states occupying characteristic output modes 
ipn (£) and <j>n Ut \z). We compute those modes and 
their occupancies (n n ) for a few realistic cases. We find 
that with increasing group velocity difference between the 
pump and Stokes waves, the scattering process involves 
an increasing number of modes which become compara- 
bly squeezed and occupied. This is also reflected in the 
photon-count statistics we compute, which changes from 
single-mode thermal statistics in the case of no group 
velocity difference, into multimode thermal statistics. 

Next, we consider reading out the atòmic polarization 
in the anti-Stokes scattering process. We show that in 
general the readout can be decomposed into múltiple in- 
dependent beam-splitter like transformations. Each of 
them mixes the statistics contained in a pair of char- 
acteristic input field (í) and atòmic (z) modes, 
producing two characteristic output modes field $^ ouí 5 (í) 

and atòmic &n Ut \z). This fact stems from the preserva- 
tion of the total number of excitations during the anti- 
Stokes scattering. 

However, when the anti-Stokes scattering is applied to 
reading out an atòmic polarization in a particular mode, 
this consideration can be much simplified. We show un- 
der what conditions the readout process is successful in 
translating the atòmic statistics into the field. It turns 
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out that efficient readout is possible for wide range of 
group-velocity differences, however usually the pumping 
intensity must be precisely controlled. Also it turns out 
that once the readout is successful, it translates a given 
atòmic mode into a predetermined field mode, and the 
latter is unsullied by statistics of any other atòmic modes. 

Let us note that the equations governing Stokes and 
anti-Stokes scattering have closed-form anlytical Solu- 
tions Q in the limit of vanishing dispersion. In this limit 
modal decomposition can be obtained by a single-step nu- 
merical diagonalization of the analytical field and atòmic 
correlation functions 0|. 

Finally let us remark, that our treatment applies to 
the case of writinga weak quantum signal into an atòmic 
medium 



oi writmg a weak q 
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APPENDIX A: DERIVATION OF THE 
NONDEGENERATE REDUCTION THEOREM 

The Bloch-Messiah reduction has been derived for a 
general case of a monolithic quantum system plT |. The 
original theorem comprises derivation of a modal expan- 
sion of the system quantum field operator, whose appli- 
cation reduces given output-input relations into a set of 
paral·lel single-mode squeezing transformations. 

However, in the particular case of Raman scattering 
process two distinct types of operators exist — those 
labelled by space and those labelled by time. In this 
fruitful extension of this theorem can be made. A 
natural division into atòmic and light subsystems which 
are separated before and after interaction is utilized for 
this. These subsystems play symmetric roles in the pro- 
cess. This is a crucial fact we will employ to specify a 
particular form of Bloch-Messiah reduction in the non- 
degenerate case. Below we will show how to introduce 
separate modal expansions for the atòmic and light sub- 
systems and cast the input-output relations into parallel 
two-mode squeezing transformations. 

We start our generalization of the Bloch-Messiah re- 
duction theorem by first reducing the number of atòmic 
and light modes involved in the interaction to a finite 
number, so that we can use the original form of the 
theorem |2Ï|. This is easily done when one notes that 
the functions appearing in our considerations can be ap- 
proximated with any chosen accuracy on a discrete grid. 
Therefore we discretize the spatial-temporal region in 
which the interaction occurs by introducing sets of points 



in time {t n } and space {z n } and restricting ourselves to 
the rectangular grid spanned by this set. Then we can 
introduce the input and output vectors: 



àin.n = «(0, t n ) 
hn,n = b(z n , — T) 

the Green matrices: 



Óout,n — à(L, t n ) 

b ut,n = b(z n , T) 



Ca n , 



$a (jn •> %rn) 



Sb ( %n i t-rxi ) 



(Al) 



(A2) 



and the characteristic mode vectors 



^ o : t} =4 out \tn) 



'k,n 



r(out) _ ,(ouí) 



= C (A3) 



Oout 




C a 0" 








'0 S a ' 




m 


_bout_ 




C b _ 




_bin_ 


+ 


S b 







which allow for casting the equation <|HJl into a matrix 
form: 



(A4) 



c s 

where the column vectors are concatenations of their 
components, and we defined block-matrices C and S. 
Since the atòmic and light subsystems play a symmet- 
ric role in the process, the C commutes with matrix O 
defined in the following way 



O 



i 

-i 



(A5) 



while for S the following holds: OS = SO f . 

The equation l|A4Jl is a Bogoliubov transformation, 
thus the Bloch-Messiah reduction can be applied [2ÏJ to 
C and S and we can express them in decomposed form: 



C = cosh ( k 

k 

S = sinh (k 



'(out) 



~i(out) 



T(in) 

n > <t>, 



k 



(A6) 

* 

(A7) 



where [ip^ , <$L ut ^] and [ij>£ n \ 0^™'] are left and right 
eigenvectors (output and input eigenmodes), defined as 
concatenation of respective component vectors, ( n are 
real numbers. Here f denotes transposition and complex 
conjugation. By definition of singular value decomposi- 
tion [23] the singular vectors obey the following relations: 



7*(m) 



{in) 
k 



ct 



tYouí) 
0k 



= cosh^fc 
= coshCfe 



~]*(out) 



(A8) 
(A9) 



multiplying both equations by O from the left and using 



CO = OC shows, that if [iPÍ out \ $ { k out) } and [$* n) , $ n) ] 
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are left and right eigenvectors, then [iip k ,—i<jr k ] and 

[iipk , — *0fe"' ) ] also satisfy the above relations, and thus 
by definition of SVD are left and right eigenvectors of 
C. Therefore in the summation IjAfij) we can group those 
vectors together, obtaining: 




t r 



i'{in) 7{in) 



~]*(out) 



l Vk , -w k 



^ OMÍ) vi m) 





(A10) 



Which shows that in each term of the singular value de- 
composition of C the off-diagonal elements vanish. Ap- 
plying analogous transformation to the summation l!A7j) 
it can be confirmed that the same property holds for 
S. These statements are the content of our extension 
of Bloch-Messiah reduction. They allow for expressing 
the blocks of C as in Eq. l(9ajl and i|9hjl in the text, and 
blocks of S as in Eq. ijflcjl and l|9dl) . 



APPENDIX B: PROOF OF DECOMPOSITION 
INTO BEAMSPLITTER TRANSFORMATIONS 

Analogously to Appendix El we begin by discretizing 
the annihilation operators 



^in,n — ('((), t 7í J 
d>in,n — d^Z n: T) 

the Green matrices: 



^out,n ^(-^5 tn) 
^out,n — ^(^n;^) 



^c„ -m — C c (t n ,t m ) Cd n m — Cd(z n , z m ) 

Sc nJn ~ S c {trii Zm) Sd nm — Sd{z n: t m ) 

and the characteristic mode vectors: 



(Bl) 



(B2) 



k,n 



\z n ) (B3) 



Next, let us rewrite Eq. 12311 in a matrix form: 



^out 








^in 


_^out_ 




—Sd Cd 







(B4) 



" c c 


s c ~ 




\ct 






\cl 


-S* r 




' c c 


s c ~ 




—Sd 


c d _ 




si 


c d ]_ 




sl 


c d ]_ 




—Sd 


c d _ 





For the energy to be conserved in the anti-Stokes scatter- 
ing process, the underbraced matrix U must be unitary 



1 
1 

(B5) 

therefore eight equations constrain the Green matrices. 
Let us first focus on the diagonal parts of the product. 
The condition C C C\ + S c Sl = 1 means that matrices 
C C C\ and S C S\ commute, therefore C C C\ and S C S\ have 
the same eigenvectors. This in turn imposes the condi- 
tion that the left singular vectors of C c and S c are the 
same. By analyzing the remaining equations contained in 
the diagonal portion of Eq. IIB5JI we conclude that the sin- 
gular vectors are shared between Green functions. Thus 
we can specify their singular value decompositions: 



n n 

Sc=J2 7n*!° Ut) * $l m)T S d = 5 n^ Ut) * *i m)T 
n n 

(B6) 



Let us note that all the singular vàlues can be always 
made real, by transferring the phase factor into the as- 
sociated singular vectors. Inserting these Green matrices 
in the above form into Eq. l|B5l) yields, that for every n 
the matrix built from the eigenvalues 



U„ 



(B7) 



must be unitary. But any real unitary 2x2 matrix has 
the form 



(B8) 



where 6 n is a real angle. Inserting the singular vàlues in 
the above form into Eq. dEfil and replacing sin 6 n = ^Jrfil 
we get Eq. GT 





In 




cos 8 n 


sin 9 n 




(3n_ 




— sin 9 n 


cos 9 n 
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